Communication loophole in a Bell-EPR-Bohm experiment: standard no-signaling may 

not always be enough to exclude local realism. 
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Assuming perfect detection efficiency, we present an (indeterministic) model for an EPR-Bohm 
experiment which reproduces the singlet correlations, without contradicting Bell's original locality 
condition. In this model we allow the probability distribution p\ of the state A at the source 
to depend parametrically on the orientation £ of one of the measuring devices: p\(X,£). In a Bell 
experiment, no-signaling between the source and each one of the devices would seem clearly sufficient 
to rule such an influence; however, not even schemes where the choice of observables takes place 
during the on-flight time of the particles can prevent, in some situations, a model of this type from 
violating the local bounds. In particular, a random shift p\(X,£,i) — > Pa(A,^2) p\(X,£„) 
allows the model to perform a "subensemble selection" for each of the terms involved in the inequality 
(analogous to what goes on with the efficiency loophole), whenever some correlation of those p\- 
shifts with the sequence of measurement choices is allowed. That correlation does not necessarily 
imply signaling during the photon on-flight time. 

PACS numbers: 



Previous note: We have been made aware (much 
thanks to Dr. M. Hall for that) of some recent develop- 
ments on the subject fffd - ITZl /. So far to our knowledge, re- 
sults regarding our model for the singlet correlations can 
be understood as a particular case of what Dr. Hall refers 
to as "measurement dependence", and in general as a 
particular case within the broader class of models studied 
in most of those references ( though the formal approach 
is quite different). Nevertheless, the intentions of this pa- 
per being less ambitious than most of the former, we still 
think it may retain some interest as a simple, hopefully 
easy to grasp counterexample to some widespread beliefs 
about the rigidity of the frontier between quantum and 
classical behaviours. On the other hand, it may also be 
useful to attract some more attention on the subject. As 
a general remark, we are not so interested in "measure- 
ment dependence" as a mathematical resource as we are 
in the fact that we may find a feasible and intuitive physi- 
cal counterpart for that phenomena, for instance thinking 
in the action of field lines generated by the measurement 
devices on the source. Besides, our approach to indeter- 
minism (see also additional notes) and its relation with 
Clauser and Home's factorability , as little sophisticated 
as it may perhaps seem, can nevertheless be, in our opin- 
ion, clarifying. 

Following John Bell 0, the results of the two distant 
measurements A, B £ {±1} performed in a bipartite ex- 
periment (for instance, of the EPR-Bohm type [![) can 
be expressed as 

A = A(a,X), (1) 
B = B(b, A), (2) 

where a, b are the orientations of each device, respec- 



tively, and where A £ A, with a probability density func- 
tion p\(X), summarizes the state of the pair of particles 
at the source. Bell's locality condition stands, therefore, 
simply for the fact that A[a, A) does not depend on b, 
nor does B(b, A) depend on a, either. This (determin- 
istic) Bell locality condition is the main hypothesis not 
only behind the original Bell inequality l2j[, but also (to 
our knowledge) for any other inequality |25j |. 

However, a more realistic view demands a generaliza- 
tion of the former expressions to the so-called indeter- 
ministic case 12611: 



A = A(a,X,u a ), 
B = B(b, A, u>b), 



(3) 
(4) 
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where uj a e SJ^Wi e 17 b are other hidden variables, 
with density functions p u (ui a , a), p u (wb, b), representing 
the state of the measuring device (and therefore para- 
metrically dependent on the corresponding orientation 
a, b). Now, for any function fj(A, B) (subindex standing 
for "joint") we have 

(f j {A,B)){a,b)= [ dXduj a duj b 

P(j)(X, Lo a ,u b ) ■ fj(A(a,X,u a ),B(b,X,0Jb))- 

(5) 

where ptj) (A, oj a , u>b) is a joint density function; just as we 
did before, we simply need to invoke the intuitive idea of 
locality to see that we must demand that A, uj a ,ujb are sta- 
tistically independent, i.e. their joint probability density 
is factorizable: p^(\,u a ,w b ) = p\(\)p u {u> a ,a)p u {ub,b). 

This said, in experiments testing Bell inequalities of, 
for instance and again, the EPR-Bohm type [lj , a " loop- 
hole" appears (there is room for a local model repro- 
ducing the observed results) when the two devices are 
not sufficiently far apart from one another, and a signal 
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can be thus (causally) transmitted between them, in the 
time that elapses between the two almost-simultaneous 
measurements. This is also known as the "locality" or 
"communication" loophole, and it would in principle be 
excluded by the standard "no-signaling" condition (the 
time between the measurements is not enough to allow 
for any sub-luminal transmission of information between 
the two locations). In any case, a violation of a Bell in- 
equality is only meaningful if neither that no-signaling 
nor any of the other possible additional hypothesis in- 
volved is violated too: in the following, we will provide 
a model that works as a counterexample on any of the 
former (bipartite) inequalities (a local model capable of 
reproducing the singlet correlations) , without necessarily 
violating the usual no-signaling condition. This will allow 
us to show that, in the first place, standard no-signaling 
between the two measuring devices should always be sup- 
plemented with no-signaling between the source and each 
one of the devices. But, moreover, this last restriction 
may be not be enough either. 

Let us now consider that the probability density func- 
tion for A is parametrically dependent on the orientation 
of the nearest (to the source) device, for instance, let it be 
p\(X,b). Expressions ©-((4]) remain valid, and we also 
demand that the three random variables X,to ai ujh keep 
their independence (their joint density function remains 
factorizable) , so now, for an arbitrary function fj(A,B) 
we have (from now on, except otherwise stated, the space 
of integration is, naturally, A © Vt a © ilf,): 

(fj(A,B))(a,b) = J d\dLu a diu b x 

Px(X, b)p u (u; a , a)p w (w 6 , b) ■ fj(A(a, A, u a ),B(b, A, u> b )), 

(6) 

using the fact that pu\{X,u) a ,u)b) is of course still factor- 
izable, as well as 

(f s (A))(a) = JdXduj a x 

px{X,b)p LJ (u a ,a) ■ f s (A(a,X,uj a )), (7) 
(f s (B))(b) = JdXdu b x 

■p x (X,b)p u (oj b ,b) ■ f s (B{b,X,w b )). (8) 

Now, all models of this kind satisfy CH's factorability 
condition p| [27j , as it is not difficult to prove HU , which 
enables us to write 

(AB)(a,b) = J dX- Px (X,b) x 

[ P(A = +l\a, X)P(B = +l\b, X) 
+ P(A = -l\a, X)P(B = -l\b, A) 
- P(A = +l\a, X)P{B = -1|&, A) 
- P(A = -l\a, X)P(B = A) ]. 

(9) 

At this point, we should realize that for any given 
P(A — p,\a, A), \i = ±1 satisfying the axiomatic laws 



of probability [29], one can always find p^{io a ) consis- 
tent with that choice, and the same for P{B = p\b,X) 
in relation to p u (io b ) . So, therefore, fixing P{A — 
p\a, A), P(B — p\b, A), together with p\(X,b), is all we 
need to define our model. Let us now make the choice 

P{A = p\a,X) = -[l + /xcos(a-A)], (10) 

P(B = p\b,X) = ^[l-/xcos(6-A)], (11) 

with, again, p <E and let us also consider that the 

state at the source A follows a probabilistic distribution 
governed by the following density function 

P\ (A, 6) = | [ 6(X - b) + 6(X - b + tt) ] , (12) 

where the Dirac deltas introduce a strong, but well al- 
lowed by our hypothesis after all, parametrical depen- 
dence of p\ on the direction b of one of the measuring 
devices. Using the last three expressions in ([9]) is enough 
to get to what we were looking for: 

(A-B)(a,b) = -cos(a-b), (13) 

which is precisely the quantum correlation for the singlet 
state. Obviously, p\ is not rotationally invariant, but it 
will be again capable of producing (apparent) rotation- 
ally invariant correlations from the point of view of the 
experimenter [30| . Besides, this model meets the stan- 
dard no-signaling requirement between the two observers, 
but, surprisingly enough, still remains completely local. 
A first consequence of this is that we need to include 
no-signaling also between each of the pairs (observer, 
source), if we want to discriminate between local and 
non-local models. In a highly ideal experiment where 
the whole correlation spectrum is analyzed, this may be 
enough, but perhaps not in a real one, as we will show 
here. 

But first, what if the parametrical dependence is for 
instance on £ ^ b? It is a matter of algebra to see that 
for 

Px (X,0 = I[S(X-0 + S(X-^ + tt)], (14) 
we obtain 

(A ■ B)(a, b) = - cos(a - b) + sin(a - £) sin(6 - £). 

(15) 

We want to show to what extent a model of this kind 
could work in an actual experiment; for instance, let us 
consider the case where the CHSH 3] inequality is tested, 
and let a, a' and b, b' be the two pairs of alternative ori- 
entations of the devices at each side. Moreover, let us 
adopt the procedures of Aspect's more restrictive experi- 
ment [1] , working with photons (and therefore subjected 
to the efficiency loophole, that we will ignore), but still 
interesting for us (see [Hj]) because it uses (uncorrelated) 
post-selection of observables for both particles, as it is 
described, for instance, in (20| . 
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Now, we will suppose that p\(X, £) obeys (fT4f . 
but is randomly shifting from one to another within 
a set of four possible density functions given by 
p\(\,a),p\(\,a'),p\(\,b) and p\(\,b') (therefore, £ £ 
{a, 6, a', 6'}), the four possibilities with equal probabil- 
ity. These shifts are random but not completely uncorre- 
cted with the choices of observable at the devices. This 
correlation does not necessarily need signaling (between 
device and source) during photon's on-flight time, regard- 
less of whether we use observable post-selection or not: 
however, there must exist some device-source communi- 
cation, though this can present some delay as well. 

We will measure that correlation by a certain param- 
eter r. For instance, if we are going to measure (polar- 
izations, projections of spin) a, 6, then there is a proba- 
bility r that £ is one of the chosen orientations (either a 
or b): a fraction T of the events where we measure a, b 
will contribute (to the observed value of given correla- 
tion (A ■ B)(a, b)) with a value — cos(a — b). On the other 
hand, for a fraction 1 — T of those events, either £, = a' 
or £ — b' (in principle let us suppose with equal prob- 
ability), which means that other contributions appear, 
modifying the overall observed correlation. Now, using 
(fTSf we obtain that, on average over all events, 

(A ■ B)(a,b) = -rcos(a - b) 

— 1(1 — r) [cos(a — b) — sin(a — a') sin(6 — a')] 
-|(1 - T) [cos(a - b) - sin(a - b') sin(fe - b')} 
= — cos(a — b) + ^(1 — T)[ sin(a — a ) sm{b — a') 

+ sin(a-6')sin(fr-6') ]. 

(16) 

Analogous expressions for {A ■ B) (a, b'), {A ■ B) (a', b) and 
(A ■ B)(a', V) (see (H) finally lead us to 

/8m(T) = q + |(1 - r) X 
[ sin(a — b') sm(b — b') + sin(a — b) sin(6' — 6) 
+ sin(a' - b') sin(6 - b') - sin(a' - b) sm(b' - b) ] , 

(17) 

with j3 q — — cos(a— &)— cos(a— b')— cos(a' — &)+cos(a' — b') 
the well known quantum mechanical prediction. 

Now, for T = 1, shifts in p\ and the choices of observ- 
ables 4>a £ {a, a'} and 4>b & {b,b'} are completely corre- 
lated, and we obtain the quantum value; for 1 > T > A 
p\ and 4>a, 4>b only bear some correlation but the model 
is still capable of producing a value that defies the in- 
equality, and for Y = 0.5, p\ and 4>a,4>b are completely 
uncorrelated, as a result of which the model cannot vio- 
late the inequality. 

We have seen that allowing for some dependence of 
p\ on certain parameters of the experiment (a, b, a',b'), 
together with some correlation (given by T) of that para- 
metric dependence on the choice of observables, makes 
the model still capable of producing a value going (con- 
siderably) beyond the local bound. Neither that depen- 
dence, nor that also necessary correlation, are completely 
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FIG. 1: A maximum violation of the CHSH inequality by 
the quantum prediction can be achieved, for instance, with 
a = 29, b = 9, a' = 0, 6' = 39. We have represented, for 
9 e [tt.&t]: (i) fi q = f3 m (Y = 1): * (red), (ii) /3 m (T = 0.80): 
o (green), (iii) /3 m (T = 0.50): x (blue), and (iv) the value we 
would obtain for a uniform distribution p\ (non-parametric): 
P(\-unf) = \fiq'- + (yellow). This last satisfies the inequality, 
adding a factor of | to the quantum value [32] ■ 

implausible [34| , even assuming no-signaling during pho- 
ton on-flight time, unless we design an experiment where 
an exhaustive evaluation of the whole spectrum of cor- 
relations (E : a, b — > E(a,b) for all a,b 6 [0,27r]) is per- 
formed, with a random, post-selected choice of measure- 
ments for each event. 

But, from the mathematical point of view... what is re- 
ally happening when T > 0.5? The answer is simply that, 
through the shift in the parametrical dependence of p\, 
we are allowing the model to perform what is known as 
"subensemble selection" HH : in this case the ensem- 
bles are {Ai,p\ = p\(X,^i)}. Some correlation between 
the choice of observables and shifts in p\ is needed, nev- 
ertheless, to make it possible: with T > 0.5 the full en- 
semble A = UiAi is unavoidably fairly sampled. This is 
the same mechanism that lays underneath the efficiency 
loophole, for experiments based on data rejection (see for 
instance 0, [H, H ^ ) . 

Nevertheless, it is still clear that the more exhaustive 
and restrictive the conditions of the experiment, the less 
the model can get close to the quantum prediction, so 
margins in the actual observed violations still play a key 
role. Taking this into account, perhaps the situation with 
the Clauser-Horne inequality should be studied more in 
detail elsewhere: there, violations are usually close to the 
local bound (in part due to efficiency constraints). 

Some of the first precedents of this work are Q and 
PH . Inparticular, we use the same "Malus cosine law" as 
Scully Q for our probabilities (TT0|) — (fTT j) of detection (in 
Scully's work, probabilities of passage through an Stern- 
Gerlach device). A complete analogy with the quantum 
case is not achieved there, as a difference with Barut and 
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Meystre's work [TTJ] , where the bridge is indeed built, 
although it also needs (besides some sophisticated math- 
ematics) the introduction of an additional (and a bit ob- 
scure) condition (some projector acting only for one of 
the devices). We also need a sort of additional assump- 
tion to make our model's and the quantum prediction 
meet, but in our case its interpretation comes up as per- 
fectly clear: we are talking of a causal influence of one 
of the measuring devices on the source of the state. This 



influence (taking place during the particle on-flight time 
or not) can find, without much imagination, a physical, 
quite plausible counterpart: the effect of far field lines 
generated by the device (an Stern-Gerlach, a polarizer or 
whatever it is). 

The author thanks Dr. R. Risco-Delgado for his com- 
ments, and acknowledges support from the Departments 
of Applied Physics (II and III) at Sevilla Univ. 
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based on the GHZ states [T^]) or that of Braunstein and 
Caves (two parties, N observables per party, [3). 
[26] Indeterminism was firt explicitly addressed in al- 
though it remains still perfectly compatible with the 
CHSH inequality [3] and its generalizations. A neces- 
sary remark is that it is also possible to give a definition 
of deterministic local realism in terms of probabilities, 
such as P(A = fia) = f A (X,a),P(B = (j, b ) = f B (X,b) 
for [la, fib £ {±1} ( or i equivalently, P(A — fj, a \B = 
lib, A, a, b) — P(A — /x a |A, a) and the same for B). 
That definition is, however, always reducible to one of the 
kind A = A(X',a),B = B(X',b), where the hidden vari- 
able has been now redefined from A to A' = A ffl /j, and 
where {i is an extra set of variables that also describe the 
state of the source. A proof of it goes like this: let us sup- 
pose that for a given value A of the hidden variable and a 
given set M or simultaneous (compatible) measurements, 
n different sets of results TZi of those measurements can 
occur, each one with probability p;. Then, we can define 
a new random variable [i taking values fj, — [14 for each of 
the former sets of results, so, naturally, P(n — fJ,i\X) = Pi- 
If we now build a new hidden variable A' = A © we can 
obtain a completely deterministic description of the sys- 
tem, but this time over A'. 

On the other hand, a hidden variable model built on A' 
is always factorizable in the sense of Clauser and Home 
5], while the one defined over A may certainly not be 
so. Proof of this last fact is trivial, given that for A', the 
probabilities involved - P(A\a, A'), P(B\b, A') - are now 
either unity or zero. 

Therefore, the use of probabilities to define local real- 
ism does not necessarily imply indeterminism. This in- 
determinism can, on the other hand, be understood in 
the sense of some "loss of information in the measure- 
ment process", and described with the introduction of 
new variables describing the state of the apparatus and 
any other factor that may locally affect the result of 
the measurement. If we name this new (set of) variables 
uj a ,ijJb, then we will finally have A — A(X' , a,uj a ), B — 
B(X',b,uj b ). 

[27] The model that we propose satisfies the factorability con- 
dition assumed in Clauser-Horne's inequality [5J: P(A = 
fi a , B = Hb\a, b, A) = P(A = fj, a \a, A) ■ P(B = /j,b\b, A), for 
any /j, — ±1. Though there has been some controversy 
on this hypothesis [q. fl9|, we will consider the problem 
is solved already: see [2q ]. 

[28] We begin our proof by choosing, in eq. (J7J), f s (x) = 1 for 
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x — p and f s (x) = if x 7^ p, we can obtain 
P(A = M |a) = 

dAda; a • Pa(A, ^p^^.s), (18) 



A ( a , A , lo a ) — 



P(A = fiU A|a) 
Pa(A,6)- 



Aj q ■ p„(a; a ,a), (19) 



' A(a ,X ,uj a ) = f± 

Dividing each member by px we obtain 



P(A = fM\a,X) = / doj a -p w (uj a ,a), (20) 

•/ A(a, A, u} a ) — fi 

that of course can be complemented with 



P(B = n\b,X) = f 

J E 



duj b ■ puj(ui b , b). (21) 

Now, if we also choose, for eq. ((6|, fj(x,y) = 1 for x — 
p a ,y = Pb and fj(x,y) — in any other case, we can 
obtain the joint probabilities: 

P(A — ji a ,B — p b \a, b) = 

dXdiOa dujb ■ px ( A , b) p^ (ui a , a) p^ (ut , b) , 

(22) 

and 

P(A = p a ,B = p b U A|a,6) = 

P\(X, b) ■ / du a du) b ■ p bJ (u)a,a)p UJ (u) b ,b), 

(23) 

which finally means that (A — p a ,B — p b \a,b,X) = 
P(A = (i a \a,X) ■P(B = ii b \b,X). 

[29] We need to guarantee < P(A = fi\a, A) < 1, P(A = 
+l\a,X) + P(A = -l|a,A) = 1. 

[30] When we assume (112f) , the rotational invariance is appar- 
ent because the paramenter £ in p>(A,£) follows (£ = b) 
the orientation of one of the devices throughout its rota- 
tion. 

[31] Using (|15[1 we also obtain 
(A ■ B)(a, b') 

— — cos(a — 6 ) + 5(1 — r)[ sin(a — a') sin(6' — a) 

+ sin(a-6)sin(6' -b) ], (24) 

(A-B)(a',b) 

= — cos(a' — 6) + |(1 — T)[ sin(a' — a) sin(& — a) 

+ sin(a'-6')sin(6-6') ], (25) 

(A-B)(a',b') 



■ cos(o' — b ) + |(1 — T)[ sin(a' — a) s'm(b' — a) 



[32] With (AB)(a,b) = -\ cos (a - £>) for p A uniform, it is 
trivial that we get a P(\- un f) = \Pq- 

[33] This model was initially thought for an experiment with 
massive particles, therefore devoid of the efficiency loop- 
hole. Besides, we were thinking in a EPR-Bohm scheme. 
Some recent experiments with atomic qubits follow more 
complicated schemes [2l], [2^], but, although the scope 
of our idea can perhaps be broadened to include those 
other schemes, we will leave it aside for the moment, or 
for elsewhere. What we would rather do is illustrate, in 
some general situation, how standard procedures do not 
rule out models like the one we have given. 

[34] To add some plausibility to our model, we can say that, 
let Tf ught be the photon on-flight time from the source 
to the device, and r s hift the mean time between shifts of 
observable at one of the devices, either one of the follow- 
ing two conditions are, for instance, sufficient to make 
T > 0.5 possible: (a) T shi ft > Tfu ght , (b) the probability 
of a shift in <f>A, <j>B in an interval Tfu g ht is less than i. 

[35] By "change of ensemble" or "subensemble selection" we 
refer to the fact that a Bell inequality is usually only valid 
when all the evaluated correlations (or probabilities) are 
significant of the same ensemble (or set) of particles (or 
states), i.e., they are fairly sampled on the total ensem- 
ble A. Local models can exploit the lack of detection effi- 
ciency to introduce unfair sampling (efficiency loophole) , 
and now we have seen that they can also do the same by 
exploting a random shift in the parametrical dependence 
of px. 

[36] J. DeGorre, S. Laplante, J. Roland. Phys. 
Rev. A 72, 062314 (2005). See also 
http://lanl.arxiv.org/abs/quant-ph/0507120, A model 
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for the singlet correlations equivalent to ours is given in 
Sec. II C (Theorem 4), with, for instace, px = px(X\b). 
M. Hall. Phys.Rev.Lett. 1 05, 250404 (2010). 
See also http://lanl.arxiv.org/abs/1007.5518. While our 
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